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AN ACCURATE ANALYTICAL MODEL 
WITH A REDUCED PARAMETER SET 
FOR SHORT-CHANNEL MOS TRANSISTORS 


Andrei SEVCENCO', Gheorghe BREZEANU” 


Rezumat. Este. prezentat un model analitic compact, bazat pe numai cativa 
parametri, pentru caracteristicile electrice ale tranzistorului MOS submicronic 
pentru afi utilizat in proiectarea analoga. Efecte secundare, ca viteza de saturare si 
modularea lungimii canalului sunt incluse atat pentru dispozitivele nMOS cdt si 
pentru cele pMOS. Sunt propusi noi parametri pentru descrierea vitezei efectului de 
saturatie. Este demonstrata o aproximare liniara a dependentei vitezei purtatorilor 
de sarcina de campul electric. Este realizata o buna concordanta intre model si 
datele experimentale pentru caracteristicile de transfer si de iesire si pentru 
transconductan{a, mdsurate pe tranzistoare nMOS si pMOS standard (0,18 um si 
0,35 yum). 


Abstract. A compact analytical model, based only on a few parameters, for the 
electrical characteristics of submicron MOS transistors to be used in analogue 
design is presented. Secondary effects like velocity saturation and channel length 
modulation are included in the model for both nMOS and pMOS devices. New 
parameters are proposed for the description of the velocity saturation effect. A 
linear approach of the dependence of the carrier velocity on the electric field is 
proved. A very good agreement between model and experimental data is achieved 
for transfer, output characteristics and transconductance measured on standard 
0.18 and 0.354um nMOS and pMOS transistors. 
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1. Introduction 

Technology advances have brought the microelectronic industry to never before 
seen levels of efficiency in production of ICs. This is a feat supported mainly by 
the extreme scaling of the most basic cell of electronic devices, the MOS 
transistor. This has not always been done in the best way needed. Since customer 
ease of use has always been a primordial condition for good business, in the 
beginning of the microelectronic industry as we know it today MOS devices were 
scaled through the constant voltage technique. This was due to TTL compatibility 
needs of the, then current, IC market. This in turn created very high fields in the 
channel of the MOS transistor, as channel lengths were ever smaller while 
voltages applied to these channels were always constant. 
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After reaching a point where secondary effects caused by these very high 
electrical fields could no longer be tolerated and because it was clear that bipolar 
technologies would never win against MOS transistors, the constant field scaling 
technique took hold. This brought about a real benefit for the IC business since 
now the scaling, and, consequently, the efficiency of production could be taken to 
ever increasing heights. 


After years of carefree business making.and profiting from the seemingly 
“unending” capability of silicon to absorb our needs, it has become clearer year by 
year that this has brought us to a turning point. Constant field scaling has taken 
supply voltages so low that we are nearing the limits of our science. 


While this has been true and ominous for quite some time in the digital world, 
analogue applications have gotten a stay of execution since what we need from 
MOS transistors used in most analogue devices is amplification. And since many 
analogue applications do not call for extremely complicated designs, design for 
wafer space efficiency has not been such a problem that scaling had to be taken to 
extremes. If we add the fact that the dynamic parameters of MOS transistors 
cannot be precisely guaranteed when considering using MOSFETs in highly 
advanced technologies, i.e. 65nm or so, we find that the analogue domain is still 
using channel sizes which digital devices have forgotten a decade ago. 


Still, we have gotten to a point were even in the analogue domain the feature sizes 
used cause unwanted secondary effects. Short channel sizes bring about short 
channel effects like velocity saturation or channel-length modulation, which are 
exacerbated to a very high degree. This work focuses on an analytical MOSFET 
model, including these two short channel effects. 


While striving for as accurate a description of MOS operation as possible, the 
comprehensive model presented in this paper takes into account the carrier 
saturation velocity for both carrier types, holes and electrons, and includes the 
Early effect. 


A linear approximation of the quadratic expression of the electron velocity versus 
the electrical field in the channel is, for the first time, introduced in order to obtain 
an analytical solution for the nMOS electrical characterisitics. 


Analytical equations are derived for the drain current dependence on gate and 
drain voltage and for the transconductance in strong inversion. 


An intuitive physical model of MOS electrical behaviour which uses a small 
number of parameters is obtained. 


The model is confirmed by comparison to the experimental data, measured on 
both nMOS and pMOS transistors obtained in standard (0.18, 0.35 and 1 um) 
CMOS processes. 
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2. Velocity saturation 


In analogue applications the most important short-channel effect in MOS 
transistors is, undoubtedly, the carrier velocity saturation. 


This effect stems from the increase of the longitudinal electrical field present in 
the channel. 


Scaling is now generally done at constant voltages and the electrical field in the 
channel has reached values. much greater than the critical electrical field, €.. Due 
to this fact, for submicron channels, the model with constant mobility is no longer 
valid. 


Mobility is a measure of the velocity acquired by the carrier as a result of the 
electrical field. At low longitudinal electrical fields, € (when the drain voltage, 
Vps> 18 small and/or L is large), the carrier velocity, v , increases linearly with the 
increase of the field and the slope is the surface mobility, wz. 


At high electrical fields, this slope decreases with the increase of the field, and the 
velocity reaches a maximum or saturation value, called the scattering-limited 
velocity, v 


scl * 


A first order analytical approximation of the dependence of velocity as a function 
of electrical field is given by[1]-[10] 


HS 
(14(E/8.)")" . 


V= 


where &.is the critical field, defined as &.=v,,/u. € has a value of 1V/um 
(1kV/cm) for electrons and 3V/um (3kV/cm) for holes in silicon [5][7]-[10]. 


Hence, in an nMOS device with L,, =0.5um a voltage drop of only 0.5V along 


the channel is needed to produce an average electrical field equal to the critical 
field. 


For holes, at 300K, q@=J as verified in [1] and @=2 for nMOS devices [3]. This 
difference is visible in Fig.l, where the carrier velocity versus longitudinal 
electrical field is plotted for both MOS device types. 


Eq.(1) is used in order to obtain the drain current expression of MOS transistors. 
In obtaining the solution for the drain current, the @ parameter has always been 
considered with a value of 1 for all carrier types, even though it was 
experimentally proved that a@=2 for electrons. In this paper we will, for the first 
time, take into account the differences determined by the values of a@ between 
electrons and holes. 
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Fig. 1. Carrier velocity dependence on longitudinal electrical field. 


3. The model 

The velocity saturation effect was included up until now in our other papers [7]- 
[10]. Models for pMOS (a=/) operation in weak and strong inversion were 
achieved. 


A set of formulas for the drain current, Jp, and transconductance, g,,, based on a 
very reduced set of material parameters were obtained [7]-[10]. This made the 
model extremely easy to handle on paper, what one might call a “back-of-the- 
envelope” quality. This simple model is extended for nMOS transistors in the 
present paper. 


To accommodate velocity saturation begin with the general equation for the drain 
current in strong inversion [1] 

p= wr|a| (2) 
where Q is the induced charge in the channel and W is the channel width. At a 
distance y along the channel, the voltage with respect to the source is V.s and the 
gate-to-channel voltage in this point is Vos-V.s. We assume that this voltage 
exceeds the threshold voltage V7. Thus, the induced electron charge per unit area 
in the channel is [1][2] 


|Q|(v) =C,.(Ves —Vr —Vis (yy) (3) 
where C;,, is the specific oxide capacitance. The electrical field in the channel is 


_ Wes 
ai (4) 
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where dV-,s is the incremental voltage drop along the length of channel dy at a 
distance y from the source. Substituting eqs.(1), (3) and (4) in eq.(2) gives 


Ws 


ae. (5) 
dV. s 
(ef) | 


Separating the variables and integrating along the whole length of the channel 
gives 


Ip =WIC,,.Ves —Vr —V,5(y))] 


Ox 


Ko 
ju aioe fk } | ag Wut Cox War Ves (y)]dV.s (6) 


where Vo), =Vcs —V; is overdrive voltage. 
The integral from the first term has an analytical solution for a@ = 1 (pMOS only). 


For electrons (i.e. @=2), we propose to linearize the carrier velocity expression, 
since for the pMOS case the drain current integral proved solvable. 


Therefore, we equate 


Mg Hg 


v= = SS 


(Seaver I 


where a is a linearization parameter. Solving with respect to this parameter we 
obtain 


_ ee, fp, Eoey” 
Vos.m Ln (8) 


The longitudinal electric field along the channel has been approximated by its 
average value 


S= Vos, m ( 9) 


Vpsm is the extracted value of Vps for which the error in using the linear 
expression to calculate the velocity is minimum. 


This observation is confirmed by plotting the velocity versus electric field in 
Fig.2. A good fitting is obtained between the curves based on the quadratic 
expression and the linearization equation, respectively, for a 0.13um nMOS 
device. 
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Fig. 2. Comparison of the real dependence on electrical field to linearized version of electron 


velocity. 


In practice the accuracy of the linearization is important only in the electric field 
range which corresponds to the longitudinal field values from actual submicron 
MOS devices. While focusing on practical biasing, we observed an important 
increase in linearization performance for shorter channel lengths. This was 
obvious in the drop in linearization error from a maximum of 9% in the 0.354m 
process to a negligible 3% for the 0./4am_ technology. This increase in accuracy 
can be observed in Fig.3. 
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Fig.3. Resulting error in the linearization process, for various MOS technologies. 


By using the general linearized expression of carrier velocity, we obtained an 
integrable form of eq.(6): 


L Yps ' Vos , c 
I, Tydy =(, Wk [Vov —V.s ]dV.s -{, De (10) 
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Carrying out these integrations we obtain 
= 1 Wy Vos 
Ip = n Va Fa VonVon Ye . (11) 


4. Electrical characteristics 


Eq.(11) is valid in strong inversion for the triode operating region. In order to 
determine an expression for the drain current in the active (or saturation) 
operating region, let Vos ac: be the value of Vps at the boundary between triode and 


: : . dl 
active regions of the transistor that sets —? =0. From (11) results 
DS 


Wisc _ AY 
,act 
! jae 2(Ves We) +1 oe 
2 DS,c 


where V,,. . 


for Vps in (11) we obtain the active region expression for Jp: 


=aLé, corresponds to the critical drain voltage. By substituting (12) 


==—kV 
2 L n-’ DS,act (13) 


Thus, we observe that the resulting general forms of the drain current in strong 
inversion we obtained here are almost identical to the expressions presented in the 
pMOS case in our previous papers [7]-[10]. 


The difference is contained in the multiplication of Vps,. by the linearization 
factor, a. 


Eq.(13) predicts that the drain current is independent of Vps in the active region, 
because channel- length modulation was not included. But the Early effect 
becomes important in submicron transistors. 


Consider the effect of changes in Vps on the channel length. Since drain voltage 
increases over Vps,ac, the pinch-off region of the channel is extended and this 
causes a decrease in the effective channel length, Lg. A Taylor series can be used 
for this dependence: 


~ (Vos SV iia ) 
* Gy, (Ves —Vpsat)t% LE) 1 2 ie 
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-1 
dL 
Ve =| - -L 
7 #) (15) 


is called the Early voltage of the transistor, which is a constant (i.e. independent 
of I D): 


To find the drain current in the active region affected by the Early effect substitute 
Leg given by (14) for L in eq.(13). The result is: 


a) oS DS,act 3] L DS ,act V;, 


cL _| Way 0 IW yy [1+ ve 

(16) 
The above expression evinces for the first time that the channel modulation effect 
is due to the effective voltage (Vps - Vps.act) which biases the pinch-off region. In 
classical SPICE modelling, the Early effect has always been attributed to the 
whole drain voltage, Vps [1]-[3]. Even though this was done in order to assure 
absolute convergence of simulated equations, it cannot be ignored that it can be a 
supplemental source of errors. 


5. Transconductance 


The transconductance, g,,, is the measure of the transistor’s amplifier function, 
because it reflects the transfer efficiency from input to output. Transconductance 
in short channel devices is modified significantly by velocity saturation. Therefore 
an accurate analytical description of MOS transconductance must include the 
velocity saturation effect. 


In strong inversion, transconductance is obtained from eqs.(16) and (12) by 
differentiating: 


dl D W k Vos act p tes W 4 Vos act 


Sin > Pn ee = 
dV; 1+ Dog Vas DW 2Val oe fl+ Worle.) hy) 


Rearranging with respect to the drain current, the above expression becomes 


OMhing, 
on Mp (18) 


where 
Eine = Wk ag, (19) 


is the saturated transconductance, caused by the velocity saturation effect. Ip,x, the 
knee drain current, is a new parameter of the model. The expression of this current 
results from eqs.(17) and (18). 
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2 2 (20) 


Figs.4 plot the transconductance versus [p and Voy , respectively. At low bias, g,, 
is inversely proportional to the square root of the current or overdrive voltage, 
respectively, and the velocity saturation is not significant. When the velocity is 
saturated, the transconductance reaches its maximum value given by eq.(19). gine 
is a constant independent of Jp, Voy and L. This.corresponds to high longitudinal 
electrical fields (i.e. € >€.). 
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Fig. 4. MOS transconductance relative to maximum transconductance versus: (a) drain current; (b) 
overdrive voltage 


Each of the curves from Figs.4 has two asymptotes, which correspond to the limit 
cases where the velocity saturation effect is insignificant and dominant, 
respectively. Extrapolation of these straight lines gives a flexing (knee) point 


where J, =I5, and Voy . = Sali. Note that the Vps value which corresponds to 
the knee current is still the critical drain voltage (V,,, =Vps...)- 

Hence, at high biases, over the knee point, the velocity saturation effect becomes 
predominant. 

The transconductance is almost unchanged by the channel-length modulation due 


to V, Ll V,. (see eq.(17)). 


DS ,act 


One more dynamic parameter which must be added to the small signal model is 
the output conductance. It is derived from eq. (16) 
dy Ip 


8a ayy (21) 


Zds 18 a result of the channel-shortening effect. 
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6. Theory-experiment comparison 


The validation of the model presented in the previous section is accomplished by 
comparison to experimental data measured on pMOS and nMOS transistors produced in 
0.18, 0.35 and 1 um commercial CMOS processes, which represent the classical 
technologies for analogue circuits. While these minimum channel sizes cannot be 
considered thoroughly deep submicron, it is still true that phenomena that exist in 
devices with feature sizes smaller than these values are governed by quantum 
mechanics, rather than by the classical theory. In this case, MOS device 
capabilities for analogue applications are diminished. However, even in 
submicron CMOS processes the MOS transistors in a given analogue circuit are 
deliberately designed to have larger channels for high performances. 


Figs.5 show a comparison between the model and experimental transfer 
characteristics measured on a pMOS device realised in the 0.18um technology and 
a nMOS device fabricated in the 0.35um technology. The knee drain currents 
calculated using eq.(20) are also inserted in Figs.5. A very good agreement of the 
theoretical curves, described by eqs.(12) and (13), and the measured data, has 
been achieved. 
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Fig.5. Drain current as a function of overdrive for: (a) pMOS transistor;(c) nMOS transistor 


We can discern on these plots the two limit situations of the drain current 
variation. When the saturation velocity effect is negligible, from (12) and (13) 
results 

1 


Inp=r—k - 


7 rae 


(22) 
This square law dependence of /p on Voy is evinced in Figs 5, for I, <<Ip x. 


To examine the limit case when the velocity is completely saturated, let 
big Ol nee 
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From eqs.(12) and (13) we obtain 


1 WwW 
Tp =a 7 Yes.cVov (23) 
A linear increase of drain current versus overdrive voltage can be observed at high 
currents when the carrier velocity is saturated (Figs.5), in agreement with eq.(23). 


The measured output characteristics.are.compared with the theoretical curves, 
based on eq.(16), in Figs.6. 


Only the active region is represented, where the output characteristics are linear. 
By extrapolating this part of the characteristics, back to the Vps axis, the Early 
voltage can be deduced. 


A good theory-experiment fit was achieved for the plotted characteristics from 
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Fig.6. Drain current as a function of Vps for: (a) pMOS transistor;(c) nMOS transistor 


Table 1 contains values for extracted parameters from the fit between the model 
and the measured data, for both transfer and output characteristics. 


The threshold voltage was extracted through the classical technique, using the 
square root of the drain current. 


The current factor, k’, and the linearization parameter, a, were obtained by 
comparing the transfer characteristics of the measured MOSFETs in strong 
inversion with the theoretical curve, given by eqs.(12) and (13) (Figs.5). We can 
see (Table 1) that the current factor increases at low channel lengths with the same 
rate for both pMOS and nMOS transistors. 


The value of a is multiplied in the table by the critical electrical field. While 
PMOS devices have a constant critical field for any CMOS process, in the 
NMOSFET?’s case, &. is different from the classical value [1] and it decreases 


with shorter channel length. 
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The Early voltage was extracted from the output characteristics (Figs.6). 


TABLE 1. Parameters extracted from the measured data for all studied devices 


CMOS Process 
Parameter 

nMOS pMOs 
L(uum) 1 0.35 | 1 0.35 0.18 | 0.18 
W(um) 25 0.5 25 0.5 2 0.4 
VV) 0.623. | 0.659 | -0.7 | -0.623 } -1.42.| -0.34 
k@4V’) | 160 [210 | 45 | 60 116. | 100 
V.(V) -353 | -50.7 | 55.6 | 16.3 — =| 
a*€(V/um) | 1.74 | 1.22 | 3 3 3 3 


A complete strong inversion plot for the output characteristic of a /zm channel 
nMOS device is shown in Fig.7. 


As we have seen earlier, the active region measurements agree completely with 
the model. Still, when considering the triode region, a noticeable overestimation 
of the real current in the channel is obtained with the model. 
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Fig.7. Complete strong inversion plot for drain current as a function of Vps for an nMOS transistor 


This difference, which increases at higher gate voltage (Fig.7), arises from the 
modulation of the carrier mobility in the channel. 


Due to the high electrical fields created by the gate biasing the carrier mobility is 
additionaly reduced determining an increase of the importance velocity saturation 
effect. This fact shows the importance of the inclusion of velocity saturation in 
any analytical model for the MOSFET. 
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The transconductance curves calculated using (16) are compared with the 
experimental data in Figs.8. The same good theory — experiment agreement as for 
the transfer characteristics (Figs. 5) is obtained, for any channel width (Figs.8). 
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Fig.8. Transconductance dependence on drain current for: (a) pMOS; (b) and (c) nMOS 


6. Conclusions 


An accurate analytical model with a reduced parameter set for submicron MOS 
transistors was presented. 


The model takes into account the velocity saturation in the channel and the 
channel-length modulation. 


Simple equations with just a few parameters were obtained for the drain current, 
drain voltage and transconductance for both pMOS and nMOS transistors. 


For nMOS transistors, a linearization of the carrier velocity dependence on the 
longitudinal electrical field was necessary. 


Several new parameters were introduced in the model in order to evince the 
velocity saturation importance. 


The knee point corresponds to the bias which creates the critical eletrical field. 


For biasing over the knee point the velocity saturation effect becomes 
predominant. 


A new expression of the drain current, which takes into account the channel 
length modulation, was achieved. 


A new expression of the drain current in the active region, which takes into 
account the channel-length modulation, was obtained. 


This equation evinces that the drain current increases due to the voltage drop on 
the effective pinch-off region. 
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The accuracy of the proposed model was proven by comparing it with the 
experimental data measured on pMOS and nMOS transistors fabricated in the 
standard 0.18, 0.35 and 1 um CMOS processes. 


A very good theory-experiment agreement was obtained for the transfer 
characteristics, output characteristics and the transconductance. 


It was demonstrated that the model correctly predicts the significant effect of the 
carrier velocity saturation on several analogue parameters of the MOS transistor. 


The model, being based on a reduced set of parameters, can be easily 
implemented in an analogue ICs simulator. 
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